{
volScalarField alphaEff_0_001("alphaEff_0_001",turbulence->nu()/Pr_0_001);

fvScalarMatrix TEqn_0_001
(
fvm::ddt(T_0_001)
+ fvm::div(phi, T_0_001)
- fvm::laplacian(alphaEff_0_001, T_0_001)
+ U.component(0)*gamma_0_001     //if your flow is periodic in y or z directions, e.g. component(1) or  //component(2).
);

TEqn_0_001.solve();
}
